Lensing of 21-cm Fluctuations by Primordial Gravitational Waves 
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Weak-gravitational-lensing distortions to the intensity pattern of 21-cm radiation from the dark 
ages can be decomposed geometrically into curl and curl-free components. Lensing by primordial 
gravitational waves induces a curl component, while the contribution from lensing by density fluc- 
tuations is strongly suppressed. Angular fluctuations in the 21-cm background extend to very small 
angular scales, and measurements at different frequencies probe different shells in redshift space. 
There is thus a huge trove of information with which to reconstruct the curl component of the lens- 
ing field, allowing tensor-to-scalar ratios conceivably as small as r ~ 10 -9 — far smaller than those 
currently accessible — to be probed. 
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One of the principle aims of early-Universe cosmol- 
ogy is detection of the inflationary gravitational-wave 
(IGW) background [JJ via measurement of the curl pat- 
tern [5] that it induces in the cosmic microwave back- 
ground (CMB) polarization. Likewise, a principle aim of 
physical cosmology is measurement of the distribution of 
atomic hydrogen during the "dark ages," the epoch af- 
ter recombination and before the formation of the first 
stars and galaxies, via detection of hydrogen's 21-cm line 
[3H5]. Several experiments are poised to soon detect the 
21-cm signal from the epoch of reionization and there 
are longer-term prospects to delve into the dark ages [7] . 
In this paper, we show that angular fluctuations of the 
21-cm intensity may ultimately provide an IGW probe 
that extends to amplitudes smaller than those currently 
accessible with the CMB. 

Weak gravitational lensing of galaxies by large-scale 
density perturbations [8] was detected in 2000 [9] and is 
now a chief aim of a number of ongoing and future galaxy 
surveys. These efforts seek the lensing-induced distor- 
tions of galaxy shapes. Weak lensing of the CMB by 
density perturbations was detected recently [TT] . The ob- 
servational signatures here are lensing-induced position- 
dependent departures from statistical isotropy in the two- 
point CMB correlation functions, or equivalently, the 
four-point correlation functions induced by lensing [TO"] . 

Primordial gravitational waves can likewise lens both 
galaxies and the CMB [T2TfT3] . The most general lensing 
pattern can, like the CMB polarization, be decomposed 
into curl and curl- free parts [15j . Since density pertur- 
bations produce (to linear order in the deflection angle) 
no curl in the lensing pattern, the curl component pro- 
vides an IGW probe. The problem, however, is that the 
curl signal, even with the most optimistic assumptions 
about IGWs, is well below the noise for both current 
galaxy surveys and even for optimistic next-generation 
CMB experiments. 

Here we consider lensing of intensity fluctuations in 
the 21-cm signal from atomic hydrogen in the dark ages. 



Atomic hydrogen in the redshift range 30 < z < 200 
can absorb radiation deep in the Rayleigh- Jeans region 
of the CMB [3J. Measurement of this absorption, over 
some narrow frequency range (corresponding to a nar- 
row redshift range), over the sky thus maps the spatial 
distribution of hydrogen at that redshift. The angular 
power spectrum of these 21-cm fluctuations extends to 
multipole moments I ~ 10 7 (limited only by the bary- 
onic Jeans mass) [3J, far larger than those, I ~ 3000, to 
which the CMB power spectrum extends (beyond which 
fluctuations are suppressed by Silk damping). The sig- 
natures of gravitational lensing of these 21-cm angular 
correlations are precisely the same as those of lensing of 
the CMB temperature map — local departures from sta- 
tistical isotropy. We can therefore adopt unchanged the 
mathematical formalism for lensing of the CMB. 

Our work resembles in spirit that in Ref. [16| which ar- 
gued that the huge number of Fourier modes available in 
21-cm maps of the dark- age hydrogen distribution would 
provide considerable statistical significance in detecting 
the IGW distortion to matter fluctuations. However, 
they consider the intrinsic distortion to matter fluctu- 
ations by IGWs. On the other hand, we consider the 
distortion to the images of the matter distribution by 
lensing by IGWs. Our work is related to that of Ref. [17j . 
who considered reconstruction of the lensing field due to 
density perturbations with 21-cm fluctuations. 

The most general deflection field A can be written as 
a function of position h on the sky as |15j . 

A = V S (f>{h) + V e -x fi(n), (1) 

in terms of curl-free (Vg^) and curl ( VVx Q) components. 
The angular power spectrum for the curl field f2(n) due 
to lensing of sources at redshift z by IGWs with power 
spectrum Pr(k) is 

c2 = ij j^MQ (2) 
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FIG. 1: The power spectrum for the deflection-field curl com- 
ponent for lensing of sources at various redshifts by a scale- 
invariant spectrum of IGWs of the largest amplitude (r = 0.2) 
consistent with current measurements. We also superimpose 
noise power spectra for lensing reconstruction carried out to 
various values of Z ma x. Also shown is the noise power spec- 
trum we estimate from co-adding the signals from all possible 
redshifts, assuming an Z max = 10 6 . 
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and 770 and ij(z) are the conformal time today and at 
redshift z, respectively. Here T(w) ~ Zj\(w)/w is the 
gravitational- wave transfer function, and j n (x) are the 
spherical Bessel functions. The angular power spectra 
for the lensing of sources at several redshifts are shown in 
Fig. [I] for L < 6, the source-redshift dependence is weak 
for a scale-invariant gravitational-wave background. 

We now review how this power spectrum is measured 
following the treatment of lensing of the CMB in Ref . [H] , 
focusing on a single redshift slice first. Given a map 1(h) 
of the 21-cm intensity as a function of position h on the 
sky, the minimum-variance estimator for the spherical- 
harmonic coefficients for the curl component of lensing 
is 
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where (7™ ap = (7; + Cf is the angular power spectrum 
of the map with Ci the power spectrum of the 21-cm 
intensity and Cf the noise power spectrum, and the sums 



are only over l+l'+L =odd. We use lower-case I for CMB 
fluctuations and upper-case L for the lensing-deflection 
field. Here, 
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where A e ^/ M are estimators for odd-parity bipolar- 
spherical-harmonic coefficients [18 in terms of the 
spherical-harmonic coefficients a™ p of the 21-cm map 
and Clebsch-Gordan coefficients C\^\, _ m ,. The estima- 
tor for the power spectrum of the curl component of the 

EJfQ| 2 /(2Ul). The 
variance of fl lm under the null hypothesis is given by 



deflection field is then 
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This noise power spectrum is plotted in Fig. [T] using the 
21-cm power spectra from Ref. [3J and taking the noise 
power spectrum Cf = for I < l n 



and Cf = 00 



for 



I > Z max - We show results for several Z max which are, 
roughly speaking, the maximum value of / with which 
the 21-cm power spectrum can be measured with high 
signal-to-noise. The signal-to-noise (squared) with which 
IGWs can be detected is then 



(S/Nf=J2(L + l/2)(C?) /(a2) 4 



(7) 



Before reviewing the numerical results, it is instructive 
to consider an analytic estimate of the noise power spec- 
• To do so, we use the flat-sky approximation 
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For L -C I we approximate \L — l\ ~ I — Lcosa, where 
cos a = L ■ I, and C|£_;~] — Ci — L(cosa)(dCi/dl). If 
Ci oc l n , then 



(cr£) 2 = J J da i-L 4 sin 2 a cos 2 . 
k /(64tt). 
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The flat-sky calculation is accurate for L > 20 and 
overestimates the noise by up to 30% at smaller L. As 
shown in Fig. 2 in Ref. [3J, the 21-cm power spectrum 
extends without suppression out to I > 10 6 , and values 
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of ^max ~ 10 7 are perhaps achievable with a bit more ef- 
fort. However, given the rapid suppression of the 21-cm 
power spectrum at higher I, the return on the investment 
of noise reduction in terms of higher l max will probably 
be small above i max ~ 10 7 . 

We now approximate the f2 power spectrum (for r = 
0.2) as C2 ^ 10~ n (L/2)- 6 . Although this approxima- 
tion differs from the numerical results for different red- 
shifts z at L ~ 30, it is quite accurate for all 30 < z < 200 
for the smallest L where most of the signal arises. From 
Eq. Q , the signal-to-noise with which the gravitational- 
wave background can be detected is 

(S/N) ~ 4.5 (W10 6 ) 2 (n/2) 2 (L^/2)- 1 , (10) 

where L m - m is the minimum L that can be measured. 

There are several things to note about this result: (1) 
The signal-to-noise obtained with the adopted fiducial 
values for Z max , L, and n is significant. (2) The scaling 
of the signal-to-noise with Z max is very rapid, and greater 
than what might have been expected (oc Z m ax) naively. 
The origin of this rapid scaling is similar to that for de- 
tection of the local-model trispectrum [19] (as the signal 
we are measuring here is, strictly speaking, an intensity 
trispectrum). Thus, the sensitivity to a gravitational- 
wave background increases very rapidly as the angular 
resolution of the map is improved. (3) The sensitivity 
decreases as L m j n is increased, so good sky coverage is 
important for gravitational-wave detection. 

While a signal-to-noise of 4.5 is respectable, and could 
be improved with even larger Z max , we can go much fur- 
ther: By changing the frequency at which the 21-cm map 
is made, we look at spherical shells of atomic hydrogen 
at different redshifts. Suppose, then, that we have 21- 
cm maps at two different frequencies that correspond 
to spherical shells separated along the line of sight by 
a comoving distance SR. Those two maps are statisti- 
cally independent at the highest I (where the vast ma- 
jority of the signal-to-noise for IGW detection arises) if 
(8R/R) > Z _1 . If AR is the separation in comoving ra- 
dius corresponding to the entire frequency range covered 
by the observations (say, redshifts z ~ 30 — 200), then the 
total number of statistically independent maps that can 
be obtained is N z ~ (AR/5R) ~ l(AR/R) ~ 0.15 1. If 
so, then the signal-to-noise from all these redshift ranges 

can be added in quadrature, and the signal-to-noise then 

1 1% 

increases by a factor N z . But there may be room for 
even more improvement: If most of the lensing occurs at 
redshifts z < 30 (as is the case for the lowest L), then 
the lensing pattern is the same for all redshift shells in 
which case every redshift shell contributes coherently to 
an estimator for SIlm- In this case, (cr^) 2 is decreased by 
factor iV~ , and the signal-to- noise increased by a factor 
N z relative to the single- z estimate. Since most of the 
signal comes from the lowest L, we estimate that the 
signal-to-noise for IGW detection obtained by coadding 



redshift shells will be 

(S/N) tot ~ 6.8 x 10 5 (WW 6 ) 3 (n/2) 2 {L^/2)~ l , 

(11) 

assuming (as above) the largest currently allowed IGW 
amplitude r ~ 0.2. Put another way, the smallest tensor- 
to-scalar ratio that can be detected at the 3tr level is 

r ~ 10- 6 (L min /2) (Z max /10 6 )" 3 (n/2)- 2 . (12) 

Note that the dependence on l max is very steep, and in- 
cluding all the information to Z max = 10 7 could yield a 
detection threshold of r ~ 10~ 9 . The full-sky calculation, 
including a more realistic shape of Ci , yields a result con- 
sistent with this estimate (Fig. [T]) . 

To put this result in perspective, we note that the cur- 
rent upper bound r < 0.22 comes from WMAP measure- 
ments of temperature-polarization correlations, although 
not from B-mode null searches. The forthcoming gen- 
eration of sub-orbital B-mode experiments are target- 
ing r < 0.1, and a dedicated CMB-polarization satellite 
might then get to r ~ 10~ 2 [20] . 

Measurement of gravitational-wave amplitudes r < 
0.01 with CMB polarization will have to contend with 
the additional contribution to B-mode polarization from 
gravitational lensing (by density perturbations) of pri- 
mordial E modes [21] . The two contributions (IGW and 
lensing) to B modes can be distinguished if the lensing de- 
flection angle can be reconstructed with small-scale CMB 
fluctuations [25J H3]- This may allow values r ~ 10 -3 
to be probed, although it requires a far more sophisti- 
cated CMB experiment (with far better angular resolu- 
tion) than simple detection of B modes would require. 

Further progress in separation of lensing and IGW con- 
tributions to B modes can be obtained with 21-cm mea- 
surements [T7] of precisely the type we discuss here but of 
the curl-free lensing component (due to density perturba- 
tions) rather than the curl component from IGWs. Such 
measurements, when combined with a precise CMB po- 
larization experiment, can in principle get to IGW ampli- 
tudes comparable to those we have discussed here. Mea- 
surement of the 21-cm curl component may therefore ul- 
timately be competitive for the most sensitive probe of 
IGWs, even if a sensitive CMB-polarization experiment 
is done. Furthermore, if both 21-cm observations and a 
CMB-polarization map are available, then measurement 
of the 21-cm curl component can be used as a cross-check 
and to complement a measurement from the combination 
of B-mode polarization with 21-cm lensing subtraction. 

While we have focusscd here on the dark ages, similar 
measurements can also be performed with 21-cm fluc- 
tuations from the epoch of reionization and also with 
galaxy surveys; the critical issue will be how high Z max 
can get. While the 21-cm curl component induced by 
lensing by density perturbations at second order is too 
small to be an issue |13) , a curl component may conceiv- 
ably arise since the atomic-hydrogen distribution is not 
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perfectly Gaussian due to non-linear gravitational col- 
lapse and baryonic effects. We speculate that this curl 
component will be small for the small-!/ modes at which 
the IGW signal peaks. We also imagine that the informa- 
tion from multiple redshifts may be combined to separate 
the IGW and any bias-induced signal. 

To close, we note that the measurements we describe 
will be challenging and are very futuristic compared to 
what current and next-generation experiments will ac- 
complish. Still, 21-cm cosmology is an exciting and 
rapidly developing experimental arena, for a good num- 
ber of scientific reasons [I] , and we hope that the idea pre- 
sented here provides one additional motivation to carry 
such work forward. 

This work was supported by DoE DE-FG03-92- 
ER40701, NASA NNX10AD04G, and the Betty and Gor- 
don Moore Foundation. 



[1] V. A. Rubakov, M. V. Sazhin and A. V. Veryaskin, 
Phys. Lett. B 115, 189 (1982); R. Fabbri and M. d. Pol- 
lock, Phys. Lett. B 125, 445 (1983); L. F. Abbott 
and M. B. Wise, Nucl. Phys. B 244, 541 (1984); 
A. A. Starobinsky, Sov. Astron. Lett. 11, 133 (1985). 
[2] M. Kamionkowski, A. Kosowsky and A. Stebbins, 
Phys. Rev. D 55, 7368 (1997) |arXi v : astro-ph/ 96 1 1 1 25] ; 
M. Kamionkowski, A. Kosowsky and A. Stebbins, Phys. 
Rev. Lett. 78, 2058 (1997) [arXiv:astro-ph/9609132] ; 
M. Zal darriaga and U. Seljak, Phys. Rev. D 55, 1830 

U. Seljak and M. Zal- 



(1997) arXiv:astro-ph/9609170 



darriaga, Phys. Rev. Lett. 78, 2054 (1997) arXiv:astro- 
Cabella 



ph/96091 69]; P. 
arXiv:astro-ph/0403392 



and M. Kamionkowski, 



[3] A. Loeb and M. Zaldarriaga, Phys. Rev. Lett. 92, 211301 
(2004) [arXiv:astro-ph/0312134] . 

[4] S. Furlanet to, S. P. Oh and F. Briggs, Phys. Rept. 433, 
181 (2006) [arXiv:astro-ph/0608032] . 

[5] M. F. Morales and J. S. B. Wyithe, Annu. Rev. Astro. As- 
trophys. 48, 127 (2010) [arXiv:0910.30lo| [astro-ph.CO]] . 

[6] J. D. Bowman and A. E. E. Rogers, Nature 468, 
796 (2010); G. Paciga et al (GMRT-EoR Collab- 
oration), Mon. Not. R. Astron. Soc. 413, 1174 
(2011); www.lofar.org; www.mwatelescope.org; 
web.phys.cmu.edu/ past; www.phys.unm.edu/ lwa; 
astro.berkeley.edu/ dbacker/eor; 
www .haystack. mit . edu/ast/ arrays/Edges. 

[7] S. Jester and H. F alcke, New Astron. Rev. 53, 1 (2009) 
|arXiv:0902.0493"l [astro-ph. CO]]. 

[8] R. D. Blandford et al, Mon. Not. Roy. Astron. Soc. 
251, 600 (1991); J. Miralda-Escude, Astrophys. J. 380, 
1 (1991); N. Kaiser, Astrophys. J. 388, 272 (1992); M. 
Bartelmann and P. Schneider, Astron. Astrophys. 259, 
413 (1992); M. Bartelmann and P. Schneider, Phys. Rept. 
340,291 (2001) astro-ph/9912508 . 

[9] D. J. Bacon, A. R. Refregier and R. S. Ellis, Mon. Not . 
Roy. Astron. Soc. 318, 625 (2000) astro-ph/0003008 ; 



N. Kaiser, G. Wilson and G. A. Luppino, Astrophys. J. 
556, 601 (2001) astro-ph/0003338 ; D. M. Wittman et 
al, Nature 405, 143 (2000) [astro-ph/0003014| ; L. van 
Waerbeke et al, Astron. Astrophys. 358, 30 (2000) 
|astro-ph/0002500 . 

[10] M. Zaldarriaga. U. Seljak. Phys. Rev. D59, 123507 
(1999) [astro-ph/9810257| ; U. Seljak and M. Zal- 
darriaga, P hys. Rev. Lett. 82, 2636 (1999) [astro-| 
ph/ 9810092 ; W. H u, Phys. Rev. D64, 083005 (2001) 
|astro-ph/0105117| ; Astrophys. J. 557, L79 (2001) [astro-| 
ph/0105424|; W. Hu and T. O kamoto, Astrophys. J. 574, 
566 (2002) |astro-ph/0111606j ; M. H. Kesden, A. Cooray 
and M. Kamionkow ski, Phys. Rev. D67, 123507 (2003) 
|astro- ph / 0302536 ; A. Lewis and A. Challinor, Phys. 
Rept. 429, 1 (2006) [astro-ph/0601594] ; S. Dodelson, 
F. Sch midt and A. Valli notto, Phys. Rev. D 78, 043508 
(2008) larXiv:0806.0331| [astro-ph]]. 

[11] K. M. Smith, O. Zldiii and O. Pore, Phys. Rev. D 76, 
043510 (2007) |arXiv:0705.3980] [astro-p h]]; C. M. Hirata 
et al, Phys. Rev. D 78, 043520 (2008) [arXiv:0801.0o44] 
[astro-ph]]; S. Das et al. (ACT Collaboration), Phys. 
Rev. Lett. 107, 021301 (2011) [arXiv:1103.2l24| [astro- 
ph.CO]]. 

[12] N. Kaiser a nd A. H. Jaffe, Astr ophys. J. 484, 
545 (1997) |arXiv:astro-ph/9609043] ; S. Dodelson, 
E. Rozo and A. Stebbins, Phys. Rev. Lett. 91, 
021301 (2003) [arXiv:astro-ph/0301177] ; C. Li and 
A. Cooray, P hys. Rev. D 74, 023521 (2006) [arXiv:astro^] 
ph/0604179 ; D. Sarkar et al, Phys. Rev. D77, 103515 
(2008) |arXiv:0803.1490| [astro- ph]]; S. Dodelson , Phys. 
Rev. D 82, 023522 (2010) [arXiv:1001.5012] [astro- 
ph.CO]]; T. Namikawa, D. Yamauchi and A. Taruya, 
|arXiv: 11 10.1718] [astro-ph.CO] . 

A. Cooray, M. Kamionkowski and R. R. Caldwell, Phys. 
Rev. D 71, 123527 (2005) [arXiv:astro-ph/0503002] . 
L. G. Book, M. Kamionkowski and T. Souradeep, Phys. 
Rev. D, in p ress, |arXiv:1109.2910| [astro-ph.CO], 
A. Stebbins, astro-ph/9609149 ; M. Kamionkowski et al, 
Mon. Not. R oy. Astron. Soc. 301, 1064 (1998) jastro- 
ph/9712030 . 

K. W. Masui and U. L. P en, Phys. Rev. Lett. 105, 161302 
(2010) [arXiv:1006.4T8i"1 [astro-ph.CO]]. 
K. Sigurdson and A. Cooray, Phys. Rev. Lett. 95, 211303 

(2005) [arXiv:astro-ph/0502549] ; O. Zahn and M. Zal- 
darriaga, Astrophys. J. 653, 922 (2006) arXiv:astro- 
ph/0511547] , 

A. Hajian and T. Souradeep, Astrophys. J. 597, L5 
(2003) [astro-ph/0308001] . 

N. Ko go and E. Komatsu, Ph ys. Rev. D 73, 083007 

(2006) [arXiv:astro-ph /0602099] . 

J. Bock et al. [EPIC Collaboration], |arXiv:0906.1188| 
[astro-ph.CO]. 

M. Zal darriaga and U. Se ljak, Phys. Rev. D58, 023003 
(1998) [astro-ph/9803150j ; A. Lewis, A. Challin or and 
N. Turok, Phys. Rev. D65, 023505 (2002) |astro-| 
ph/0106536. 

M. Kesden, A. Cooray and M. Kamionkowski, Phys. Rev. 
Lett. 89, 011304 (2002) [arXiv:astro-ph /0202434] . 
L. Kno x and Y. S. Song, Phys . Rev. Lett. 89, 011303 
(2002) [arXiv:astro-ph/0202286 . 



[13 
[14 
[15 



[17 



r 

[18 

[19 
[20 
[21 



